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ABSTRACT 

In this paper we study the colour-magnitude relation (CMR) for a sample of 172 
morphologically-classified elliptical and SO cluster galaxies from the ESO Distant Clus- 
ter Survey (EDisCS) at 0.4 < z < 0.8. The intrinsic colour scatter about the CMR is 
very small ((omt) = 0.076) in rest-frame U — V. However, there is a small minority of 
faint early-type galaxies (7%) that are significantly bluer than the CMR. We observe 
no significant dependence of a- m t with redshift or cluster velocity dispersion. Because 
our sample is strictly morphologically-selected, this implies that by the time cluster 
elliptical and SO galaxies achieve their morphology, the vast majority have already 
joined the red sequence. The only exception seems to be the very small fraction of 
faint blue early- types. Assuming that the intrinsic colour scatter is due to differences 
in stellar population ages, we estimate the galaxy formation redshift zp of each clus- 
ter and find that z-p does not depend on the cluster velocity dispersion. However, zp 
increases weakly with cluster redshift within the EDisCS sample. This trend becomes 
very clear when higher redshift clusters from the literature are included. This suggests 
that, at any given redshift, in order to have a population of fully- formed ellipticals and 
SOs they needed to have formed most of their stars ~ 2-4 Gyr prior to observation. 
That does not mean that all early-type galaxies in all clusters formed at these high 
redshifts. It means that the ones we see already having early-type morphologies also 
have reasonably-old stellar populations. This is partly a manifestation of the "progen- 
itor bias", but also a consequence of the fact that the vast majority of the early-type 
galaxies in clusters (in particular the massive galaxies) were already red (i.e., already 
had old stellar populations) by the time they achieved their morphology. Elliptical and 
SO galaxies exhibit very similar colour scatter, implying similar stellar population ages. 
The scarcity of blue SOs indicates that, if they are the descendants of spirals whose 
star-formation has ceased, the parent galaxies were already red when they became SOs. 
This suggests the red spirals found preferentially in dense environments could be the 
progenitors of these SOs. We also find that fainter early- type galaxies finished forming 
their stars later (i.e., have smaller zp), consistent with the cluster red sequence being- 
built over time and the brightest galaxies reaching the red sequence earlier than fainter 
ones. Combining the CMR scatter analysis with the observed evolution in the CMR 
zero point we find that the early-type cluster galaxy population must have had their 
star formation truncated/stopped over an extended period At > 1 Gyr. 
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1 INTRODUCTION 

Galaxy clusters have proven to be very useful laboratories 
for the study of galaxy formation and evolution. They can 
provide large and diverse galaxy samples across practically 
small areas of sky. Although the relative importance of na- 
ture and nurture in shaping galaxy evolution remains de- 
bated, it is well established that many galaxy properties in 
the nearby Univers e correlate strongly with their environ- 
ment. In particular, IPostman fc Gellerl |l98J) and iDressler] 
(1980) found that the morphological mix changes as a func- 
tion of local galaxy density, with early-type galaxies (E/SO) 
dominating in locally-dense environments, while late types 
(spirals and irregular galaxies) become increasingly impor- 
tant towards lower densities. Moreover, at low redshift there 
is an obvious bimodality in the colours of cluster galaxies 
th at is well correla t ed with morpholog y (e.g. see recent work 
bv lConselice]|2006l ; IWang et al. 1120071 ). The "blue cloud" is 
dominated by spirals and irregular galaxies, and the promi- 
nent ridge of red galaxies (the so-called "red sequence" ) is 
mainly populated by E/SO galaxies. This implies that mor- 
phology correlates at some level with the stellar population 
content (and hence colour). 

iBauml J 19591 ). iFabed j 19731 ) and ICaldwelll (|l983l ) 
established the existence of a red sequence of clus- 
ter ellipticals in the local Universe and showed that 
these galaxies have s ystematically redder colours wi th in- 
creasing luminosity. Visvanathan fc Sandagei l| 19771 ) and 
ISandage fc Visvanathanl 1 1978al lbl) found that this colour- 
magnitude relation (hereafter CMR) is universal. The de- 
tailed study of the UV K colours of local clust e r early -type 
galaxies carried out by iBower. Lucev fc Ellis! (J1992) con- 
firmed Sandage & Visvanathan's anticipation that both SOs 
and ellipticals follow the same relation. Furthermore, they 
also showed that the observed scatter about the CMR is 
very small (~ 0.04 mag in U — V for their sample). 

In the past decade, a number of studies have shown 
that the CMR of elliptical galaxies holds at p rogressively 
higher rcd shifts, at least up to z = 1.4 (e.g. [Ellis et al. 



1997; Stanford et al. 1998; van Dokkum et al. 2000, 200 



Bla keslee et al. I 120031 ; iMei et al. I 120061 ; iLidman et al 



20081 ) . As a consequence, the CMR is, arguably, one of the 



most powerful scaling relations obeyed by the early-type 
galaxy population at the cores of clusters, encoding impor- 
tant information about their formation history. 

The slope of the CMR has traditionally been interpreted 
as the direct consequence of a mass-mctallicity relation (e.g. 
iFaberl Il973l ; iLarsonl [l97i ; iGallazzi et "al~l 120061 ) . The clas- 
sical explanation of this mass-metallicity sequence is based 
on the idea that star-formation-induced galactic outflows 
would be more efficient at expelling metal-e nriched gas in 
low-mass g alaxies t han in massive ones ( e .g. 'Lar son 1974; 
Dekel fc Silkl 1 19861: iTremonti et al~l [2004 : Ee Lucia et al. 
2004a; KobavaihT Springe! fc Whitell20071; iFinlator fc Dave 



2008l : lArimoto fc Yoshiilll987l ). An alternative interpretation 



in which the CMR is predominantly an age sequence would 
imply that the relation changes significantly with redshift 
as less massive galaxies approach their formation epochs. 
The possibility that age is the main driver for the CMR was 
ruled out by observations of clusters at intermediate redshift 
that sho wed that the slope of the CMR evolves little wit h 
redshift (JKodama fc Arimotol Il997l ; iKodama et al. 1 1 19981 ). 



Nevertheless, weak age trends along the CMR have been 
claimed (e.g . Ferreras et al. 1 1 19991 ; iPoggianti et al. I l200ll ; 
iNelan et al. 1 120051 ) , even though it seems clear that they 
are not the main physical driver . 

IBower. Lucev fc Ellisl (|l992T ) interpreted the small scat- 
ter about the CMR as the result of small age differences at a 
given galaxy mass. The tightness of the relation then implies 
very synchronized star-formation histories for these galax- 
ies. Larger colour scatter would imply later episodes of star- 
formation, or a wider range in galaxy formation redshifts. 
These results are not only indicative of the passive evolution 
of elliptical ga laxies but also of an early formation epoch 
(z > 2-3; e.g. IBower Lucev fc Ellisl Il992l ; [Blakeslee et al.l 
l2003uMei et al. 1120091 ). Studies of absorption-line indices in 
the spectra of earl y-type galaxies also imply old ages (e.g. 
iTrager et al. Ill998l ). Some evidence has been found that the 
mean stellar ages of early-type galaxies may depend on their 
stellar mass in the sense that lower-mass galaxies appear to 
have formed t heir stars at later ep ochs than t he more mas- 
sive o nes (e.g. iThomas et al. II2005I ). although ITrager et al. I 
( 2008) find no such trend in their study of the Coma cluster. 

The above interpretation of the nature of the CMR, al- 
though traditionally accepted, has an important problem: 
it assumes that all red-sequence galaxies that we see today 
can be identified as red-sequ ence members of high r edshif t 
galaxy clusters. As noted by Ivan Dokkum fc Franxl |l996j), 
this assumption is probably wrong because of the so-called 
progenitor bias: if the progenitors of some early- type galaxies 
were spirals at higher redshift, they would not be included in 
the higher redshift samples, which biases the studied popu- 
lation towards older ages. This effect has b een corroborated 
by rec ent studies of the CMR evolution. |Pe Lucia et al. I 
|2007j) found a significant deficit of faint red cluster galax- 
ies at 0.4 < z < 0.8 compared to galaxy clusters in the 
local Universe. They conclude that the red sequence po- 
pulation of high redshift clusters does not contain all the 
progenitors of nearby red sequence cluster ga laxies (see also 
|Pe Lucia et al. II2009I . and references therein) . iTanaka et al. I 
(1998) also find such deficit in a galaxy cluster at Z — 1.1. 
We will come back to this issue in Sections [4] and [7] 

In this paper, we study the CMR for a total sample of 
174 morphologically-selected elliptical and SO galaxies con- 
tained in 13 galaxy clusters and group s at 0.4 < z < 0.8 
from the ESO Distant Cluster Survey ([White et al. 1120051 . 
EDisCS), taking advantage of the availability of extensive 
Hubble Space Telescope (HST) imaging obtained with the 
Advanced Camera for Surveys (ACS) and extensive ground- 
based imaging and spectroscopy (see Section [5}. We inter- 
pret the scatter about the CMR as a proxy for age (or 
formation time ii?), and study its dependence on intrinsic 
galaxy properties such as their luminosities and morpholo- 
gies, and the role of the environment as quantified by the 
mass/velocity dispersion of the clusters. We complement the 
scatter analysis with information derived from the zeropoint 
of each cluster's CMR to constrain not only the formation 
epoch of early-type galaxies but also the duration of their 
formation phase. 

It is important to point out that, even though there is 
much evidence suggesting that the CMR scatter is princi- 
pally driven by galaxy age ( e.g. Kodama fc Arimotq Il997l ; 
iKauffmann fc Charlotl 1 19981 ; iBernardi et al. 1120051). metal- 



licity variations could also contribute to it 1 Nelan et al. I 
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Figure 1. HST F814W images for a representative sample of ellipticals (top row), SOs (second row) and "blue tail" galaxies (third row). 
The last row shows the 2 galaxies that were excluded from our study due to morphological misclassification. The rightmost 2 galaxies in 
the third row (blue tail) exhibit some degree of perturbation in their morphologies (see text for further discussion). The horizontal white 
lines correspond to 5kpc. 



120051 ). If that is the case, the stellar ages we derive in this 
paper would be lower limits since the colour scatter we mea- 
sure would contain both an ag e component a nd a metallic- 
ity one. Taking the results of iNelan et al. I (|2005l ) at face 
value, the maximum scatter in metallicity at a given ve- 
locity dispersion (luminosity) is ~ O.ldex, implying that at 
most we could be systematically underestimating the stel- 
lar ag es by ~ 0.15dex, where we have assumed IWorthevI 
II 19941 ) 3/2 age-metallicity degeneracy law. Notwithstanding 
this caveat, even if the absolute ages were affected at this 
level, it is not unreasonable to expect that the effect on rel- 
ative ages (the main focus of this paper) would be smaller. 

This paper is organized as follows: our dataset and the 
sample of early-type galaxies are described in Section [5] In 
Section [3] we explain the "scatter-age method" used in this 
work to study the formation histories of early-type galaxies. 
Section U shows the scatter measured about the CMRs, Sec- 
tion [5] the results of the scatter-age test and Section [6] the 
zeropoint analysis. Finally, we discuss our results and draw 
conclusions in Section [7] 

Throughout this work, we use Vega magnitudes and 



adopt a WMAP cosmology l|Spergel et al. 1120071 ) with Qm = 
0.28, fi A = 0.72 and H = 70Kms" 1 Mpc" 1 . 



2 THE DATA 

EDisCS originated as an ESO Large Programme designed to 
study cluster structure and c luster galaxy evolu tion over a 
large fraction of cosmic time ([White et al. 1120051 ) . The com- 
plete dataset consists of 20 fields containing galaxy clusters 
at redshifts between 0.4 and 1. The sample was constructed 
selecting 30 of the highest surface brightne ss candidates in 
the L as Campanas Distant Cluster Survey (JGonzalez et al. I 
120011 ). and confirming the presence of an apparent cluster 
and a possible red sequence with VLT 20-min exposures in 
two filters. From these candidates, 10 of the highest surface 
brightness clusters were followed up in each of two bins at es- 
timated redshifts 0.45 < z CBt < 0.55 and 0.75 < z cst < 0.85, 
where z cst was based on the magnitude of the putative 
brightest cluster galaxy. 

For the 20 fields with confirmed cluster candidates, 
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Table 1. Properties of the cluster sample. The columns correspond to: cluster ID, spectroscopic redshift, line-of-sight velocity dispersion, 
number (N) of early-type galaxies (E+SO), number ellipticals (iV(E)) taken into account in the scatter calculation, and numb er of 
"blue tail" members (- / V^iue) - See text for details. Cluster redshifts and velocity dispersions were taken from lHallidav et al. I (120041 ') and 
iMilvane-Jensen et al.l J2008I) . 



Cluster name 
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^(km/s) JV(E+S0) AT(E) N hlue 



C11037.9-1243a 

C11138.2-1133a 

C11138.2-1133 

C11232.5-1250 

C11354.2-1230a 

C11103.7-1245a 

C11227.9-1138 

C11054.4-1146 

C11040.7-1155 

C11054.7-1245a 

C11054.7-1245 

C11354.2-1230 

C11216.8-1201 



0.4252 
0.4548 
0.4796 
0.5414 
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deep optical imaging was taken using the FORS2 camera 
and spectrograph at the Very Large Telescope (VLT). A 
full description of the available optic al images and photom- 
etry is given in IWhite et al. I (|2005l ). In brief, the optical 
photometry consists of B, V and / imaging for the 10 in- 
termediate redshift (z cst — 0.5) cluster candidates and V, 
R and / imaging for the 10 high redshift (z cs t ~ 0.75) 
cluster candidates. Typically, the integration times were 
45min for the intermediate redshift sample and 2h for 
the high redshift sample. In addition, near-IR J and K 
photometry was obtained using the SOFI camera at the 
New Technology Telescope (NTT) (Aragon-Salamanca et 
al., in preparation). This was obtained for all but three 
fields, two of them were rejected due to bad weather 
and an additional field was rejected because its spectro- 
scopic redshift histogram was not cluster-like. FORS2 at 
the VLT was also used to obtain multi-obj ect spectroscopy 
for ~ 100 objects in each field (sec Halliday et al. 1 [20041 ; 
iMilvang- Jensen et al. 1120081 ). The redshift data were used 
to confirm the physical reality of the clusters and measure 
their velocity dispersions. This process revealed that several 
of the fields co ntained m ultiple clusters and gro ups at dif- 
feren t redshifts ( Halliday et al. 1200 4 ; Milvang- Je nsen et al. I 
|200S| ; iPoggianti et al. 1120061). Additional fol low-up includes 
narrow-band Hcv imaging (iFinn et al. 120051 ') and XMM X- 
ray observations (|Johnson et al. 1 12006) for a subset of the 
clusters. Moreover, ten of the highest redshift clusters in the 
sample were also imaged with HST's ACS camera in the 
F814W filter, from which reliable visual morphologies were 
derived for a larg e sample of cluster and field galaxies (see 
iDesai et al. 1120071 . for details). 

For consistency with our previous work on the CMR of 
EDisCS cluster galaxies (|De Lucia et al. M2007I '). in the fol- 
lowing we use magnitudes and colours measured in seeing- 
matched images with FW HM = 0.8 arcsec (the typical 
seeing in the optical images), using a fixed 1.0 arcsec ra- 
dius circular aperture. This aperture represents a compro- 
mise between minimizing sky-subtraction and contamina- 



tion errors and being as close as possible to measuring global 
colours. We note that 1.0 arcsec corresponds to ~ 7kpc at 
z ~ 0.7, and at these redshifts early-type galaxies in the lu- 
minosity range considered here have half-light radii ~ 3 kpc 
l|Treu et al. ll2005l ; lTfuiillo et al. 1120071 ). The associated pho- 
tometric errors were derived by placing empty apertures in 
regions of the image without detected objects to estimate 
accurately the sky contribution to the error budget. This is 
justified since the sky noise is the dominant source of error 
when measuring the apertu re magnitudes of our faint galax- 
ies (see lWhite et al. 1120051 . for details), 

In this paper, we focus on a sub-sample of the EDisCS 
dataset consisting of 174 galaxies in 13 clusters and groups. 
This selection was based on the following criteria: 

(i) The galaxies must b e spectroscopically- 
confirmed cluster/group m embers (jHallidav et al. I 120041 : 
IMilvang- Jensen et al. I 120081 ). This ensures a very clean 
sample , avoiding the uncertainties introduced by other 
cluster membership criteria such as ph otometric redshifts 
JRudnick et al. 1120091 : [Pello et al. 1120091 ). The penalty is, of 
course, a significant reduction in the sample size. 

(ii) They must have early-type morpholog y (E or SO), 
based on visual classification from HST images IjDesai et al. I 
l2007f ). We note that by imposing this, the sample reduces 
to galaxies within the 10 fields observed with HST, 

(iii) The galaxies should belong to clusters/groups with 
at least 4 early-type members in order to measure the CMR 
scatter with reasonable accuracy. This allows us to detect 
the presence of an intrinsic colour scatter with > 3u confi- 
dence in all cases. 

Since all spectroscopically-confirmed members in the HST- 
covered area have been morphologically classified, the selec- 
tion function of our sample is the same as that of the spec- 
troscopic sample. This means that for all practical purposes 
our early-type galaxy sample beha ves like the original I- 
band selected spectroscopic sample (JMilvang- Jensen et al. I 
120081 ). and therefore as a rest-frame B-band selected sam- 
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pie. On average we typically reach Mb < —18.5, with some 
cluster-to-cluster variation. Although all our analysis has 
been carried out using this empirically-defined J-band se- 
lected sample, thus maximizing the sample size, we have 
checked that using a rest-frame B-band luminosity selected 
sample would not have altered any of our conclusions. We 
have also checked that the spectroscopic selection function 
does not affect our conclusions. We produced Monte Carlo 
realizations of the colour-magnitude diagram of our sample 
taking in to account the empirical selec tion function deter- 
mined by iMilvang- Jensen et al. I (|2008l) and found that the 
simulated colour scatter agrees very well with the measured 
one. 

To test the robustness of the morphologies for the galax- 
ies in our sample, we re-examined visually their HST images 
(see Figure [1} and found onl y two galax ies that had been 
misclassified as early- types in lDesai et al. I 1J2007I ). The first 
one, EDCSNJ1138096-1135223, is clearly not an elliptical 
and shows a very perturbed morphology. The second one, 
EDCSNJ1138127-1134190, is in a dense group of (proba- 
bly) interacting galaxies and there is a brighter elliptical 
very close to the position of this object. It is obvious that 
the wrong galaxy was classified. The last two lines of Ta- 
ble [4] show some of the properties of these two galaxies and 
their HST F814W images are shown in the bottom row of 
Figure Q] It is not surprising that both galaxies are substan- 
tially bluer than the red-sequence. They are also quite faint, 
where visual classification is, perforce, less reliable. These 
two misclassified galaxies were removed from our sample and 
will not be discussed further. We noticed that 6 out of the 
remaining 172 galaxies (3.5%) have signs of perturbation, 
although their early-type morphology is clear. Interestingly, 
2 out of these 6 slightly perturbed galaxies are significantly 
bluer than the CMR. We will come back to this in Section[4] 

The cluster sample with the corresponding redshifts, 
line-of-sight velocity dispersions (a v ) and number of early- 
type members i s show n in T able [l] The a v 's are tak e n from 
lHallidav et al~l \2004 ) and IMilvang- Jensen et al. I (|2008h . 
The reliability of these velocity dispersions as mass estima - 
tors has been confirmed by weak lensing (jClowe et al. Il200q ) 
and X-ray (jJohnson et al. 1120061 ) estimates. 



3 METHOD: THE SCATTER-AGE TEST 

T he scatter-age test carried out in this paper was developed 
bv lBower. Lucev fc Ellis! |l992j) as a reasonably simple, yet 
powerful method for constraining the formation history of 
early-type galaxies. They applie d it at z ~ to galaxies 
in the Virgo and Coma clusters. Ellis et al. I (| 199 71 ) applied 
the same test to galaxies in three z ~ 0.54 clusters. We 
apply it here to a much larger cluster sample, covering a 
significant cluster mass range. We have the added advan- 
tage that since 1997 the uncertainty in the cosmological 
parameters, and thus the transformation of redshift into 
look-back time, has decreased considerably. We also ben- 
efit from a large and homogeneous galaxy sample in a wide 
range of redshift and cluster velocity dispersion (or clus- 
ter mass). This method has been used by many authors 
in th e past fe.g. |van Dokku m et al Ill9 98: Stan ford et al. I 
1 19981 ; iBlakeslee et al. I l2003t 120061 ; iMei et al. I l200d 120091 
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Figure 2. The R—I vs. I colour-magnitude diagram for the early- 
type galaxies in cluster CL1216.8— 1201 at z = 0.79 is shown here 
as an example. Elliptical galaxies are represented by "+" signs, 
and SOs by open diamonds. The solid line shows a linear fit to 
th e colour-magnitude rel ation with the slope value determined 
bv lDe Lucia et al. I 1120071 ). See text for details. The dotted lines 
correspond to ±0.3 mag from the CMR. For reference, the typical 
sizes of the error bars are plotted on the top of the figure as a 
function of magnitude. 



Irlilton et al. 1120091 ). A description of the specific steps we 
took to perform the scatter-age test follows. 

We first constructed colour-magnitude diagrams 
(CMDs) of the early-type galaxies in each cluster using 
the photometric bands closest to rest-frame U and V. 
Rest-frame U — V measures the strength of the 4000A 
break and is therefore a very age-sensitive broadband 
colour (see Section 13.11 for a detailed justification of our 
choice of observed colour). For most of the clusters this 
choice required CMDs of R — I versus /, but for the 
three lowest redshift clusters we used V — I versus I (see 
Table [2}. We then fitted a linear CMR for each cluster 
using a fixed slope of — 0.09 [j and setting the zeropoint 
to the median colour. This procedure is very robust, in 
particular for groups and clusters with small numbers of 
galaxies where the CMR slope cannot be determined to 
sufficient accuracy. Our results are not sensitive to the 
exact choice of slope since in general th e CMR is rea sonably 
flat and redshift-independent (|Holden et al. 1120041) . As an 
example, the CMD for the early-type galaxies in cluster 
CL1216.8-1201 is shown in Figure [5] 

For each cluster, the observed scatter in the galaxy 
colours about the CMR (<T b s ) was computed as the r.m.s. of 
the residuals (in the colour direction) between the observed 
colours and the fitted CMR. We reject outliers in this pro- 
cess by imposing the condition that galaxy colours should be 
within ±0.3 mag from the CMR. While other methods such 
as the biwe ight scatter estimator used by other authors (e.g. 
IMei et al. i 2009) reject outliers implicitly, we chose to do it 
explicitly. We discuss these outliers in some detail later. 



1 This value was previously used in |Pe Lucia et al. I 1120071 ) to 
construct V-I CMRs for colour selected red-sequence galaxies 
from the EDisCS database 
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Figure 3. Time-line of the Universe illustrating the different pa- 
rameters used in Equations[T]and [2] The time arrow starts on the 
left at the beginning of the universe (z = oo). The orange region 
shows the total available time galaxies can use to form stars. Ai is 
the time galaxies actually spend forming stars. From then on, the 
cluster galaxies are assumed to evolve passively until the observed 
redshift (z c [). We define £p as the time elapsed from the epoch 
when star formation ended until the cosmic time corresponding 
to z cl . 



Following iBower. Lucev fc Ellis! |l992), the intrinsic 
scatter ((Tint) was then obtained by subtracting, in quadra- 
ture, the mean value of the photometric colour uncertainty 
from the observed s catter. The colour uncertainties range 
from 0.012 to 0.021 (| White et al. 1120051 1. and have therefore 
littl e effect on the observed scat ter (see Section 2]) • 

IBower. Lucev fc Ellia |l992j) used (Tint to constrain the 
formation history of the galaxies by assuming the following 
relationship between the colou r scatter and the colour evol u- 
tion of the stellar population (jBower. Lucev fe Ellislll992h : 

8{U - V) = ^ 7 V) ° (fa - t F )ft < (T in t, (1) 

at 
where tn is the age of the universe at the cluster redshift, £f 
is the look-back time from then to the epoch at which star 
formation ended, and d(U—V)o/dt (where the subscript "0" 
denotes rest-frame) is derived from galaxy evolution mod- 
els. This factor should be reasonably well understood as it 
is mainly governed by main-sequence evolution (for a given 
IMF). In this equation, ft parameterizes the spread in for- 
mation time Ai as a fraction of the total available time: 

Thus, ft = 1 implies no synchronization, i.e. the galaxies in 
the sample formed using all the available time, while ft — 0.1 
would mean a high degree of synchronization, with all the 
galaxies forming in the last 10% of the available time. Figure 
[3] illustrates the time- line defined by Equations Q] and [2] 

We calculated the d(U — V)o/dt factor using 
iBruzual fc Chariot! J20051 . hereafter BC03) modelfl for a 
passively-evolving stellar populations that formed in a sin- 
gle burst of 0.1 Gyr duration. The exact burst duration 
does not have a significant effect on our conclusions pro- 
vided that it is much s horte r than ip. Solar metallicity 
(Z Q = 0.02), a lChabrierl l|200&t ) initial mass function (IMF) 
and no du st attenuation were assumed. Using alternative 
IMFs (e.g. ISalpeterl ll955T ) would not alter our conclusions 



2 For the range of ages discussed in this paper and the optical 
colours on which we base our conclusions, using alternative popu- 
lation synthesis models such as those of .Maraston, (2005 ) would 
not change our results. 




F (Gys) 



Figure 4. Illustration of the colour-scatter method and the as- 
sociated random (a) and systematic (b) uncertainties for the 
cluster CL1216. 8-1201 at z = 0.79. The rate of colour change 
(dColour / dt) as a function of £p is shown by the identical solid 
red lines in panels (a) and (b). In this example we use the observed 
(R — I), very close to rest-frame (U — V) at the cluster redshift, 
computed using BC03 stellar population models for galaxies with 
solar metallicity and a single star formation burst of 0.1 Gyr dura- 
tion (see text for details). The scatter about the CMR provides an 
upper limit to the allowed rate of colour evolution, parametrized 
by Equation[3]for a given ft (Equation[2}- This constraint is shown 
in both panels by the identical solid black lines, as derived from 
the intrinsic colour scatter for this cluster and assuming three 
different values of ft. The intersection between the observational 
lines (solid black) and the model ones (solid red) constrains tp. 
The dotted black lines in panel (a) correspond to the ±1<t random 
errors affecting the solid black line as a result of the observational 
uncertainty in the colour scatter. The red dotted lines in panel (b) 
illustrate the effects of systematic model uncertainties (e.g. metal- 
licity) on ip . These lines correspond to models with the same star 
formation history as for the red solid line but different metallici- 
ties: the upper line has ^sub-solar = 0.008, whilst the lower line 



has Z, 



super — solar 



0.05. 



because for the stellar masses of interest (given the range 
in £f ), the relative differences in the IMFs are only minor 
(cf. IBower. Lucev fc Ellis 1992]). The use of models with 
solar met allicity is partially motiv a ted b y the results pre- 
sented in ISanchez-Blazquez et al. I (|2009T ) , where ages and 
metallicities are derived for EDis CS red-sequence galaxies 
from their absorption line indices. ISanchez-Blazquez et al. I 
(2009) found solar-metallicity models agreed well with the 
observed galaxy spectra. We discuss the effect of assuming 
different metallicities later. 

Figure [4] shows how Equation [l] can be used to con- 
strain the star-formation history of the galaxies from the 
colour scatter o"i n t of our richest cluster. It also illustrates 
the effect of the relevant random and systematic uncertain- 
ties. In Figure [4f a), the red solid line represents d(R — I)/dt 
(very close to rest-frame d(U — V)o/dt) as a function of £f, 
calculated from the BC03 models. The black solid lines are 
defined by the equation 

d(R-I) _ 



(Tint 



(3) 



dt (t H - M/3 

for several values of ft. Equation [T] implies that the allowed 
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region lies below the black lines, and thus the intersection of 
the red and black lines provides a constraint (upper limit) 
on £f for a given /3. It is clear that the colour scatter alone 
cannot be used to constrain £f and fi simultaneously. The 
dotted black lines in Figure Ufa) represent the ±lcr ran- 
dom error^j in the colour scatter, showing their effect on the 
£f uncertainty. The effect of systematic uncertainties, such 
as changing the model chemical composition, are shown on 
panel (b). The dotted red lines correspond to stellar popu- 
lation models with non-solar metallicity (Zsub-soiar = 0.008 
and ^ S upcr-soiar = 0.05). It is immediately apparent that 
the effect of systematic model uncertainties such as these 
in the calculation of £f is, in general, significantly larger 
than that of photometric random errors. This implies that 
absolute values of tv must be interpreted with great cau- 
tion. However, it is not unreasonable to assume that these 
systematics would affect all the galaxies similarly, making 
any differential or comparative study precise and, hopefully, 
robust. Unless otherwise stated, we consider random uncer- 
tainties only when discussing £f since our study is largely 
comparative, but it is important to bear in mind that sub- 
stantial systematic uncertainties do exist. 



3.1 Colour dependence of the derived £f 

Previous studies have noted that colours which bracket 
the 4000A break provide the most sensitive indicators of 
age changes, yet the least affected by photometric errors 



CL 1216.8-1201 



(e.g. Blakeslee et al. Il200q). For this reason, and following 
iBower. Lucev fe Ellis! ( 19921 ) . we decided to carry out the 
age-scatter test using observed colours close to rest-frame 
U—V. Nevertheless, it is instructive to study how the actual 
colour choice could affect our results. We used our richest 
cluster, CLI2I6.8— 120f, as the ideal test bed for this pur- 
pose. Using the galaxy sample in this cluster, we carried out 
the scatter-age test with CMDs compiled for different sets of 
colours. Most colours straddle the 4000A break at z = 0.79 
(the cluster redshift) , with the exceptions of I — J and I — K. 
The values of <Ti n t and £f derived for each colour are plotted 
in Figure \5\ for different values of /3. It is clear that, at least 
for j3 ^ 0.3, all the colours provide a consistent value of £f 
within the combined random and systematic errors. How- 
ever, for 13 — 0.1 the optical-optical and optical-infrared 
colours give discrepant results. This is probably because at 
the relevant stellar population ages (1-2 Gyr) asymptotic 
giant branch stars have a potentially large, and very uncer- 
tain, contribution to the near-infrared galax y emission , mak- 
ing the model predict ions very unreliable IjMarastonl 120051 ; 
IConrov fc Gunnl 120091 ). The colour with the smallest scat- 
ter and smallest scatter uncertainty is R — I, which is the 
closest to rest-frame U — V at this redshift. This provides 
additional justification for the use of observed colours that 
are the closest match to rest-frame U — V in our analysis. 



3 The errors in the colour scatter were estimated from the 16% 
and 84% confidence levels in the \ 2 distribution of the measured 
Ci n t, which correspond to ±lcr uncertainties. 




0.05 



0.10 0.15 

<5 lnt (colour) 



0.20 



Figure 5. The derived formation times of the early-type galaxies 
in the cluster CL1216.9— 1201 at z = 0.79 are plotted against the 
intrinsic CMR scatter measured for the different colours used in 
the scatter analysis. With the exception of I — J and I — K all 
colours straddle the 4000 A break at the cluster redshift. Filled 
stars correspond to /3 = 1.0, diamonds to /3 = 0.3 and circles to 
/3 = 0.1 (see Equation \2$ . The black dotted line indicates the age 
of the Universe (in) at the cluster redshift. In this Figure the 
errors in ip are calculated as the sum (in quadrature) of the ran- 
dom and systematic uncertainties discussed in the text. Within 
the errors, most colours yield consistent values of tp (particu- 
larly for /3 ^ 0.3). However, for /3 = 0.1 the optical-optical and 
optical-infrared colours give discrepant results possibly due to the 
difficulty in modeling the contribution of asymptotic giant branch 
stars. 



4 THE SCATTER IN THE 

COLOUR-MAGNITUDE RELATION 

4.1 The CMR scatter for different clusters 

For each cluster, we calculated the intrinsic colour scatter 
<7int following the method described in Section [3] The values 
of (Tint and the actual colours used for each cluster are listed 
in Table[2] At this stage, we exclude from this calculation all 
galaxies whose colours are > 0.3 mag bluer than the fitted 
CMR. The number of galaxies excluded in each cluster is 
listed in Table [l] In total, there are 12 (7%) of these blue 
early-type galaxies in our sample. We justify this approach 
and discuss these galaxies later. 

The colour scatter shows no significant evolution with 
redshift for the clusters in the EDisCS sample (Figure [5] 
top panel) . To extend the redshift baseline and compare our 
results with previous studies, we plot in Figure [6] similar 
colour scatter measurements from the sources listed in Ta- 
ble[3] No redshift dependence is found even for this extended 
redshift range. The bottom panel of Figure [6] further shows 
that the colour scatter does not correlate with cluster veloc- 
ity dispersion (<r„) either, implying that the scatter is not 
strongly affected by cluster mass. We note that the veloc- 
ity dispersion range spanned by our sample is very broad 
(200 < o v < 1200 km/s). Adding the clusters in the com- 
parison samples reinforces our result. 
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Table 2. Values of the measured intrinsic scatter (a- ln t) and the calculated formation time (4p) and redshift (zp) for solar metallicity and 
/3 = 0.1,0.3 and 1.0. The associated uncertainties correspond to the random errors (16% and 84% confidence levels in a x 2 distribution). 
Recall that ip is the lookback time, from the cluster redshift, since star formation ceased. 



Cluster name 


z 


CMR 


Cint 


t P (/3 = 0.1) 


tp(/3 = 0.3) 


t F (p = 1.0) 


z F (/3 = 0.1) 


2 F (/3 = 0.3) 


zf(P = 1-0) 








colour 




(Gyr) 


(Gyr) 


(Gyr) 








C11037.9- 


-1243a 


0.4252 


V - 


- I 


09+ 003 
u ' ua -0.01 


9 q+0.4 
z -°-0.6 


4.4±°; 3 


6.ot ; 2 3 


0.76lg;0 8 


1.3 ±0.1 


2 1+ - 1 

^•1-0. 2 


C11138.2- 


-1133a 


0.4548 


V - 


-I 


10+ 003 


2 4+°- 2 


3.34±°i 


4.oi°;i 


n o 9 +0.04 
U - 8/ -0.06 


1 Qq+0.03 

i- uc> -o.05 


1.22 ±0.03 


C11138.2- 


-1133 


0.4796 


V- 


-I 


o.o9±°:°! 


9 1+0.3 


q 7 +0.2 
°-'-0.3 


5 1+0-2 


0.791°;°? 


1 in+0.07 
i - iJ -0.09 


1.8 ±0.1 


C11232.8- 


-1201 


0.5414 


R - 


-I 


0.08±°;° 2 


1 a+03 
l-S-0.4 


q O + 0.2 

J -°-0.3 


5 1+°' 2 


0.82±°;°o? 


l.40t°;°° 8 


2.6 ±0.2 


C11354.2- 


-1230a 


0.5952 


R- 


-I 


0.04i°;° u9 


22 +0.5 


42 +0.4 


c 7+0.3 
°-'-0.7 


i.oi8i 


-i 7+O.2 


9 O+0.3 


C11103.7- 


-1245a 


0.6261 


R - 


-I 


0.08±g;S 


1 o+°- 6 


3.1+"' 4 


4 7 +o.3 

4 -'-0.8 


0.8i°i 


1 4+0-2 


n 9+O.2 
z - z -0.5 


C11227.9- 


-1138 


0.6357 


R- 


-I 


09+ 004 


5+°- 4 
U.0_ 4 


o c+0.3 
z '°-0.6 


4 1+0.2 
^•-"--0.4 


71+ 007 

u -' 1 -0.06 


1 2+ - 1 
1 - z -0.2 


1 s+o-i 


C11054.4- 


-1146 


0.6972 


R- 


-I 


06+ 002 
u - UD -o.oi 


1 0+°' 3 
i ' u -0.6 


2 9+°' 2 
z ' 3 -0.3 


4 5+0-2 
4 -°-0.3 


o.seiH 7 


1 44+0-09 


2.4 ±0.2 


C11040.7- 


-1155 


0.7043 


R - 


-I 


(15+ 002 

u - uo -o.oi 


1 l+°' 4 


q 9 + 0.3 


4 9+ - 2 


0.9 ±0.1 


1.6 ±0.2 


9 7+O.2 
z -'-0.3 


C11054.7- 


-1245a 


0.7305 


R- 


-I 


H+ 006 


5+° 5 
U.0_ 4 


1.6t°' 5 


4 0+° 2 

4 - U -0.5 


80+ 01 

u - ou -0.07 


1 4+0-1 
1 '*-0.2 


2 1+ ' 1 
z - x -0.3 


C11054.7- 


-1245 


0.7498 


R- 


-I 


in+ 002 
u - iu -o.oi 


fi+°- 3 
U.O_ .4 


1 S+ ' 3 


O fi + 0.1 


n O4+0-06 


i qn+0.08 


l.90t°;° 6 


C11354.2- 


-1230 


0.7620 


R - 


-I 


n5+ 003 
u - uo -o.oi 


n q+0-5 
u - M -0.8 


2 q+0.4 
z ' y -0.9 


4 fi+° 3 


9+ - 1 
U.M_ 2 


1 fi+o- 2 

l-°-0.4 


9 7+O.3 
■'•'-0.5 


C11216.8- 


-1201 


0.7943 


R- 


-I 


n n^s+ - 009 


9+°- 2 
u - 9 -0.3 


9 S+0.2 


4 5+0-1 
4 -°-0.2 


96+ 005 
u - bD -0.07 


1.6 ±0.1 


9 7^+0.09 
z -' d -0.2 



Table 3. Main characteristics of the comparison samples. 



Cluster name 



Colour used 
in o- int 



Ref. 
for a:„ 



(Km/s) 



Ref. 

for a v 



Symbol in 
Fig.E] 



low-z 



Coma 

Virgo 

CL1358+62 

CL10412-65 

CL10016+16 

CL10054-27 



0.0231 

0.0038 

0.3283 

0.510 

0.546 

0.563 



0.056 ±0.01 
0.044 ± 0.01 
0.079 ±0.01 
0.131 ±0.027 
0.06 ±0.01 
0.06 ±0.01 



U -V 
U -V 
B-V 



(U -V 
(U -V 
(U -V 



z=0 
z=0 
z=0 



82llf 8 
632l« 
1027±|J 

68ii 2 r 5 

230 ± 18 



* 
* 



MS 1054-0321 0.831 0.070 ± 0.008 (U - B 

RX J0152. 7-1357 0.834 0.050 ± 0.005 (U - B 

CL1604+4304 0.897 0.031 ± 0.003 (U - B 

high-z CL1604+4321 0.924 0.043 ± 0.006 (U - B 

RDCS J0910±5422 1.106 0.060 ± 0.009 (U - B 

RDCS J1252.9-2927 1.237 0.112 ±0.022 (U - B 

RX J0849+4452 1.261 0.070 ± 0.014 (U - B 

RX J0848+4453 1.270 0.049 ± 0.027 (U - B 



z=0 
z=0 
z=0 
z=0 
z=0 
z=0 
z=0 
z=0 



1156 ± 82 

1203+^3 
703 ± 110 
582 ±167 
675 ± 190 
747^ 
740^ 
650 ± 170 



-74 
-84 
1-113 
-134 



very high-z XMMU J2235. 3-2557 1.39 

XMMXCS J2215. 9-1738 1.46 



0.055 ±0.018 J-Ks 

0.12 ±0.05 z 8 50-^ 



5 762 ± 265 

6 580 ± 140 



□ 
A 



Note: Following Bla keslee et al. I 1120061). the scat ter in (U — B) z —p was tra nsfor med into (U — V) z —o s catter by adding 0.04. 
References: (lllBower. Lucev fc Ellij lll992ri: f2')lvan Dokkum et al~l lll998l): (3 ) lEllis et al~l Jl997h ; (41 iMei et al~l J2009T) and references 
therein; (5) iLidman et al. I J2008T ) ; 761 Iffilton et al. I d2009| N |; (71 iBorgani et al. I dl999T) 



4.2 CMR scatter dependence on galaxy properties 

To explore the overall behavior of the galaxy colours around 
the CMR, Figure [7] shows the distribution of the residuals 
for the complete galaxy sample as a function of the absolute 
rest-frame B magnitude Mb (fo r details of the calculation 
of Mb, see lRudnick et al. 1120091 1. By construction, the resid- 
uals are concentrated about their median value (~ 0; solid 
line in Figure 0. The vast majority of the colour residu- 



als follow a normal distribution reasonably well. The scatter 
is small, as discussed above. However, at faint magnitudes 
there is a clear "blue-tail" containing a few galaxies with sig- 
nificantly bluer colours (smaller blue symbols in Figure [7} . 
These results can also be seen in the upper histogram of Fig- 
ure [8] which shows the distribution of the colour residuals 
for the complete sample of early-type cluster galaxies. The 
measured scatter for the whole sample is small (<7 bs = 0.078 
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Figure 6. The upper panel shows the dependence of the intrinsic 
scatter in the observed colour closest to rest-frame (U — V) z =o for 
EDisCS clusters (filled black diamonds) and several comparison 
samples at lower and higher redshift (grey symbols). The low 
re dshift cluster sample ( asteris k s) was compiled from the work 



of Bower, 



lucev U Ellid jl992h : Ivan Dokkum et al. I jl99Sh and 
Elhs_et_aL] Jl997|). The higher reds hift sample was take n from 
Mei et al. I (l2009h (open diamond s), Iffilton et al. I (l2009h (open 



triangle) and lLidman et al. I 1120081 ) (open square). See Table[3]for 
details about the comparison samples. This plot reveals there is 
no significant CMR scatter evolution with redshift up to z < 1.5. 
The lower panel shows the scatter as function of cluster velocity 
dispersion. In both panels, the solid line represents the median 
o"i nt value for the EDisCS clusters and the dotted lines correspond 
to ±2cr. The CMR scatter does not correlate with cluster velocity 
dispersion. 



when excluding the "blue-tail" ) and significantly larger than 
the scatter due to the photometric errors (~ 0.017). This im- 
plies that the intrinsic colour scatter for the complete sample 
is (Tint = 0.076lo'oo4i very close to the average scatter for the 
individual clusters (((Tint} = 0.077, cf. Tabled). 

There are 12 galaxies (7% of the total sample) with 
colours > 0.3 mag bluer than the CMR. The number of 
"blue" galaxies in each cluster is listed in Table [T] while 
their individual IDs and some observed properties are pre- 
sented in Table [4] These galaxies were excluded from our 
scatter-age analysis (Section[5]) for consistency with previous 



Figure 7. Observed colour residuals around the CMR as a func- 
tion of Mg of all the galaxies in the sample. Note that we express 
the magnitude as Mg — 5 Log h, where h = Hq/70. Elliptical 
galaxies are represented with crosses and lenticulars with dia- 
monds. The solid line indicates the location of the median, dotted 
lines correspond to 2<r and 3(X for the black data points, and the 
blue symbols represent the galaxies in the "blue tail" (see text for 
details). The median luminosity is also shown for reference (black 
dashed line). 



studies (e.g. iMei et al. 1120091 ) where outliers are rejected im- 
plicitly. We prefer to exclude them explicitly but discuss the 
implications that their existence and properties have in our 
conclusion. Interestingly, 2 out of the 12 blue galaxies show 
a small degree of disruption in their morphologies, as found 
by visually inspecting their HST images (see comments in 
Tableland Figure[T}. In Section[2]we found 6 galaxies (4%) 
in our full sample that, despite their clear early-type mor- 
phology, show some signs of disruption. Now we find that 2 
of these disturbed galaxies have "blue" colours, indicating 
that among the blue galaxies, morphological disturbances 
are much more common than among the ones in the CMR. 
Nevertheless, the rest of the blue galaxies (10) do not show 
any clear sign of morphological disruption. We further dis- 
cuss the implications of these faint blue galaxies (i.e. the 
"blue tail") in our conclusions. 

In what follows, we concentrate on the remaining 160 
galaxies whose colours are within ±0.3mag from the CMR. 
From these, we constructed sub-samples according to differ- 
ent galaxy properties: 

• Morphology (E vs. SO), as indicated by different sym- 
bols in Figure [7] 

• Luminosity (Luminous vs. Faint) , divided at the median 
rest frame B absolute luminosity (corresponding to Mg e = 
— 19.8; cf. vertical dashed line in Figure [7}. 

The four bottom histograms of Figure [8] show the distri- 
bution of the colour residuals for each one of these sub- 
samples. Within the errors, we find that both ellipticals and 
SOs show the same scatter. However, the luminous galaxies 
have a slightly smaller intrinsic scatter than the faint ones. 
The values of a lri t for each sub-sample are listed in Table [S] 
In the next section, we interpret these scatters in terms of 
the star-formation history of the different galaxy samples. 
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Table 4. Properties of the "blue tail" galaxies. The "excluded" 



laxies are also listed in the two bottom lines. 



EDisCS galaxy ID 



Type Residual (mag) Mg Comments (from spectra; morphology) 



EDCSNJ1054199- 


-1146065 


E 


-0.32 


-18.97 


EDCSNJ1054207- 


-1148130 


SO 


-0.41 


-18.43 


EDCSNJ1232307- 


-1249573 


E 


-0.90 


-19.51 


EDCSNJ1232336- 


-1252103 


SO 


-0.52 


-19.24 


EDCSNJ1232304- 


-1250391 


E 


-0.42 


-20.03 


EDCSNJ1138050- 


-1132546 


E 


-0.60 


-18.17 


EDCSNJ1138068- 


-1132510 


E 


-0.68 


-18.73 


EDCSNJ1138034- 


-1133049 


E 


-0.64 


-18.07 


EDCSNJ1354073- 


-1233336 


E 


-0.45 


-18.86 


EDCSNJ1354022- 


-1234283 


E 


-0.68 


-19.15 


EDCSNJ1037564- 


-1245134 


SO 


-0.37 


-17.61 


EDCSNJ1138135- 


-1137137 


E 


-0.66 


-18.28 



EDCSNJ1138096-1135223 
EDCSNJ1138127-1134190 



-0.33 
-0.33 



Several emission lines 

Considerable [Oil] emission 

Absorption-line spectra 

Galaxy of spectral type k+a 

Considerable [Oil] emission; Some sign of disturbance 

Strong [Oil] emission 

Starburst?; High surface brightness 

Strong [Oil] emission; Some sign of disturbance 

Strong [Oil] emission 

Strong [Oil] emission; Compact galaxy 

Considerable [Oil] emission 

Considerable [Oil] emission 



-19.00 Strong [Oil] emission; Very disturbed/Merging 
-17.88 Either HII regions or large merger 



* These galaxies where misclassified as E in lDesai et al. I (2007). We excluded them from the sample since their HST images reveal that 
they are not early- type galaxies and they show strong signs of disruption (see Section [2j . 
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Figure 8. Histograms of the observed colour residuals around 
the CMR for the different galaxy sub-samples computed from 
Figure [7] (a) the complete sample (Ellipticals and SOs); (b) El- 
lipticals only (c) SOs only; (d) luminous galaxies and (e) faint 
galaxies. The dotted lines in the top panel correspond to 2a and 
3ct for a Gaussian distribution with cr b s = 0.078. The open blue 
part of each histogram corresponds to the "blue tail" (see text 
and Figure [7] for details). 



5 STAR FORMATION HISTORIES 

5.1 Star-formation histories of the early-type 
galaxies in each cluster 

Using the method discussed in Section [3] and the intrinsic 
colour scatter measured for the early-type galaxies, we com- 
puted, for three values of /3, the formation times tp and 
their respective errors for each individual cluster or group 
(see table [2| . An inspection of this table immediately shows 
that, as expected, for higher values of f3 (less synchronous 
galaxy formation) older ages are required to explain the 
small colour scatter. A clear trend is also apparent: at a 
fixed /3, higher redshift clusters have smaller tF. However, 
if we correct for the difference in look-back time using our 
adopted cosmology, these £f can be translated into forma- 
tion redshifts, zf, and the trend disappears. All the EDisCS 
clusters yield consistent formation redshifts for their early- 
type galaxy population (zf — 0.8 for f3 = 0.1, zp ~ 1.4 for 
= 0.3 and z F ~ 2.4 for /3 = 1.0). This is shown in Fig- 
ure [S] where z-p is plotted vs. cluster redshift for the differ- 
ent values on /3. For comparison, we overplot the for mation 
redshifts derived by ISanchez-Blazcmez et al. I (|2009l ) using 
absorption-line indices in the spectra of EDisCS early-type 
red-sequence galaxies (blue and red diamonds). The blue 
points (corresponding to galaxies with velocity dispersions 
< 175km/s) agree very well with our formation redshifts for 
j3 = 0.3 (black diamonds), while the red points (galaxy ve- 
locity dispersions > 175km/s) agree with j3 ^ 0.3. Moreover, 
galaxy ages derived from the analysis of the Fundamental 
Plane of these galaxies (Saglia et al. in preparation) are also 
in agreement with our formation redshifts for j3 ^ 0.3. 

Figure [9] also shows that zp may be slightly higher for 
higher redshift EDisCS clusters (fitted line). Although this 
trend is not very significant, it would be desirable to extend 
the redshift baseline to test whether it continues at higher 
redshifts. W c can do that by using the study published by 
iMei et al. I (J2009T ) . These authors follow a very similar proce- 
dure to ours, and predict formation times based on a colour- 
scatter analysis for a value of j3 ~ 0.3. Their galaxy samples 
also contain morphologically-classified ellipticals and SOs, 
and can therefore be compared to ours. Their results are 
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Figure 9. Galaxy formation redshift zp vs. cluster redshift. Sym- 
bols correspond to the different values of /3 used: triangles for 
j3 = 0.1, filled diamonds for /3 = 0.3 and stars for /3 = 1.0. The 
solid line is a linear fit to the solid black diamonds. Fo r com - 
parison, the predicted zp from Sanchez-Blaz quez et al. I (2009) 
are shown as larger coloured symbols. The blue diamonds corre- 
spond to the morphologically-selected sample of EDisCS early- 
type galaxies with galaxy velocity dispersions < 175km/s, while 
the red symbols correspond to the sample with galaxy velocity 
dispersions > 175 km/s. The blue points agree very well with 
our zp for /3 = 0.3, whilst the red points are in agreement with 
our results if /3 ^ 0.3. We observe an increase of zp with cluster 
redshift (see also Figure [TUT) . 



plotted as open diamonds in Figure 1101 together with our 
results for f3 = 0.3. If we take these points at face value, the 
trend of increasing zf with redshift becomes very significant. 
However, a word of caution is required. Although our study 
and that of lMei et al. I (j2009f ) are very similar, there are some 
differences. First, their photometry and the colours that they 
use are different because the higher redshift of their clusters. 
Second, they estimate the intrinsic colour scatter using bi- 
weight scale estimator (which implicitly excludes outliers). 
However, using their method with our data does not change 
our results since we exclude outliers explicitly, as discussed 
in section [3] Finally, their implementation of the scatter- 
age test differs in s ome minor details f rom o urs, although 
they use the same iBruzual fc Charloj IJ2003I ) models. We 
believe these differences are probably not important for this 
exercise, but we cannot be completely certain without re- 
analysing their data. With all these caveats, the trend ob- 
served in Figure [TD] would imply that that morphologically- 
classified elliptical and SO galaxies formed earlier in higher 
redshift clusters than in lower redshift ones. 

Figure [11] shows zf plotted against cluster velocity dis- 
persion. Consistent with Figure [6l no correlation is found 
between zf and a v , implying that the formation time of 
morphologically-classified ellipticals and SOs does not de- 
pend strongly on cluster mass. 



Figure 10. Galaxy formation redshift zp (derived using = 
0.3) vs. cluster redshift for E DisCS clusters (so lid diamonds) and 
for the clusters published bv lMei et al. I (120091) (open diamonds). 
We observe that zp increases slightly with cluster redshift for 
the EDisCS sample (see also fitted line in Figure |9j. This trend 
becomes stronger when we include the higher redshift comparison 
sample. 
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Figure 11. Galaxy formation redshift zp vs. cluster velocity dis- 
persion. As in Figure[9] symbols correspond to the different values 
of /3: triangles for /? = 0.1, filled diamonds for f3 = 0.3 and stars 
for /3 = 1.0. The dashed, solid and dotted lines correspond to the 
median formation redshift for each value of /3. We find that zp 
does not depend on the cluster velocity dispersion. 



5.2 Dependence of the star- formation histories on 
galaxy properties 

In table [S] we show the colour scatter and derived formation 
redshift for the galaxy samples divided in terms of morphol- 
ogy and luminosity, as discussed in Section l4~2l Within the 
errors, ellipticals and SOs show the same scatter. Since the 
E and SO samples have very similar mean redshifts, simi- 
lar colour scatters imply similar average formation redshift. 
However, the faint galaxies seem to exhibit a larger scatter 
than the bright ones. This, together with the fact that the 
luminous subsample has a higher average redshift than the 
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Table 5. Scatter analysis results for the morphology- and luminosity-split sub-samples. The colour scatter and the derived formation 
times and redshifts are given for three values of 0. 



Sample 


(*> 


Cint 


t F (0 = O.l) 


t F (0 = O.3) 


i F (/3 = 1.0) 


z F (/3 = 0.1) 


z F (/3 = 0.3) 


z F (0 = 1.0) 


All 


0.6411 


07? +0.005 
uu ' ' -0.004 


2 64+ - 07 
z '° -0.08 


4 04+0- 06 
■*- u *_0.07 


5 4^+0.06 


i 99 +0.3 
x - z -0.3 


i ok+0.04 
J-^-O.OB 


9 qo+0.07 
z -»»_0.08 


E 
SO 


0.6443 
0.6589 


(VTE+0.00S 
U - U '°-0.004 
n n 7 n+0.009 
u - u ' y -0.007 


2.65±°;° L 9 

2 6+ 01 
z -°-0.2 


4 QK+0.08 

4 0+ - 1 


5 45+0.07 
°'*°-0.09 
5 o+O.l 


i 99+O.3 
x - z -0.3 
i 9 4 +0.05 


i ocr+0.04 
I ' 85 -0.05 
i o 7 +0.06 
l-»'_0.07 


9 qo+0.07 
^■»»_0.08 

9Q+0.1 

Z - M -0.1 


Luminous 
Faint 


0.7016 
0.6071 


070+°'°° e 
uu ' u -0.005 

n nsn+o- 007 


2.636 ±0.1 
2.668+. ,;! 


4.00l°f 
4.10t°;f 


(- OO + 0.08 

°' o>5 -0.09 

5 fi +0.1 
°- D -0.1 


I Q4+0.04 
1 - >5 ^-0.0B 

i 11+0.03 
± - 1± -0.04 


9 n o+0.08 
z - uo -0.09 
i fi7 +0.05 
J " D '-0.05 


n 1+0.2 

°-0.2 

2 fi^+0° 8 
Z - DO -0.09 



faint one, suggests that the most luminous (massive) early- 
type galaxies formed earlier than the fainter (less massive) 
ones. 

An identical conclusion is rea ched if the sample is split 
by stellar mass (derived following iBell fc de Jong||200ll ) in- 
stead of by luminosity, which is not surprising since homo- 
geneous colours imply near-constant stellar mass-to- light ra- 
tios. When splitting the sample by galaxy velocity dispersion 
or dynamical mass (cf. Saglia et al. in preparation), similar 
trends are observed, albeit with smaller statistical signif- 
icance since only half of the galaxies in our sample have 
measured velocity dispersions. 



6 THE CMR ZEROPOINT 

Our CMR scatter analysis cannot constrain t F and /3 si- 
multaneously. Until now, we have not discussed the zero- 
point of the CMR because predicting absolute colours 
from stellar population models is, arguably, more uncer- 
tain than predicting differentia l colour changes (see, e.g., 
lAragon-Salamanca et al. I 1 19951 ). However, since accurate 
zero points are available, it is worthwhile checking whether 
consistent and, perhaps, additional constraints can be ob- 
tained from them. Figure [121 shows for each of the EDisCS 
clusters the CMR colour corresponding to the median ab- 
solute B magnitude of the sample, after correcting for 
luminosity evolution. Specifically, for each cluster we de- 
termine the colour of its CMR fo r M B = Mb + 1.15 , 
where M^ is em pirically-determined (jCrawford et al. 1120091 ; 
iRudnick et al. 1 12009| ). These zero-points do not correlate 
with intrinsic cluster properties such as their velocity dis- 
persions or masses. 

The solid lines in th e Figure show the predictions of 
iBruzual fc Charlotl (|2003l ) models for z F ~ 1.5 (correspond- 
ing to /3 = 0.3, cf. SectionGnj for Z sub _ S oiar = 0.008 (blue), 



Zsoiar = 0.02 (black) and Z S u P cr- S oi ar = 0.05 (red). It is 
clear that for median luminosity galaxies non-solar models 
do not provide an acceptable fit to the galaxy colours for any 
value of /?|j However, solar- met allicity models do a reason- 
able job. This provides additional justification for the use of 
solar metallicity models in our analysis (cf. section [3]). 

Taking the solar models at face value, a constraint on fi 
can be derived from Figure [121 The black-dashed line corre- 
sponds to j3 — 0.1 (z F ~ 0.9), the black solid line to ft — 0.3 



4 For clarity, we show the non-solar metallicity lines for = 0.3 
only. 
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Figure 12. V - I CMR zero-point vs. redshift for the EDisCS 
cluster galax ies (points) compa r ed wi th population synthesis 
models from IBruzual fc Charlotl (120031) . The black lines have 
been computed with Z so i ar = 0.02 metallicity, the blue one with 
^sub-solar = 0.008, and the red one with Z supcr _ solar = 0.05. 
The black dashed line corresponds to /3 = 0.1 (z F ~ 0.9), the 
black, blue and red solid lines to = 0.3 (z F ~ 1.5) and the 
black dotted line to /3 = 1.0 (z F ~ 2.5). See text for details. The 
observed points clearly rule out a = 0.1 scenario (i.e. very syn- 
chronized formation for all the galaxies). However, they are in 
agreement with ^ 0.3. 



(zf ~ 1.5) and the black dotted line to — 1.0 (z F ~ 2.5). 
The observed points clearly rule out j3 = 0.1 (i.e. very syn- 
chronized formation for all the galaxies). They are in rea- 
sonably good agreement with /3 = 1.0, but /3 = 0.3 is not 
ruled out. Hence, it seems reasonably safe to conclude that 
/3 ^ 0.3 on the basis of this analysis. If we translate this into 
the time interval At over which all the galaxies "formed", 
the constraint translates to At > 1 Gyr. Since t F refers to 
the time at which star formation ceased, this implies that 
there was an extended epoch over which cluster galaxies had 
their star formation truncated/stopped. In other words, this 
cessation of star formation was not synchronized for all the 
cluster early-type galaxies. 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the colour-magnitude relation 
(CMR) for a sample of early-type galaxies from the ESO Dis- 
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tant Cluster Survey (EDisCS). Our sample consists of 172 
strictly morphologically-classified ellipticals and SO galaxies 
in 13 clusters and groups with redshifts 0.4 < z < 0.8 and ve- 
locity dispersions 200 < cr v < 1200km/s. All these galaxies 
are spectroscopically-confirmed cluster members, and their 
magnitudes span the range —22 < Mb — 51og/i < —17.5. 
We have analyzed the colour scatter about the CMR and 
its zeropoint to derive meaningful constraints on the forma- 
tion history of these galaxies. Assuming that the intrinsic 
colour scatter about the CMR is due to differences in stellar 
population ages, our main results are: 

• In agreement with previous studies, the intrinsic colour 
scatter <Ti nt about the CMR in rest-frame U — V is small 
({ffint} = 0.076). However, there is a small minority of faint 
early-type galaxies (7%) that are significantly bluer than 
the CMR and where excluded from the scatter analysis. 
These galaxies probably represent a population of young 
galaxies that have not yet joined the red-sequence popu- 
lation. Interestingly, only 2 out of the 12 blue galaxies 
show signs of morphological disturbances and/or interac- 
tions, while the rest are bona-fide ellipticals or SOs. How- 
ever, the vast majority of the blue galaxies have emission 
lines in their spectra indicative of ongoing star formation 
(see table [4]). Faint blue low-mass ea rly-type galaxies have 
been reported in prev ious studies (e.g. lBlakeslee et al. 1120061 ; 
iBamford et al. 1120091 ). preferentially in low density environ- 
ments. To explain the existence of these low- mass blue early- 
type galaxies with normal E/S O morphologies in the field 
iHuertas-Companv et al. I (|2010l ) propose two possible expla- 
nations. First, they could be the result of minor mergers, 
which would trigger centrally-concentrated star-formation, 
helping to build a bulge, and eventually taking them to the 
red sequence. Alternatively, the disks in these galaxies are 
perhaps being (re)built from the surrounding gas, moving 
then back (or staying) in the blue-cloud. It is hard to see 
how this second possibility would work in the cluster envi- 
ronment, where it is more likely that gas is removed than 
allowed to fall onto these low mass galaxies. In clusters, mi- 
nor mergers remain a possibility, in particular if they occur 
while these galaxies were in filaments and/or groups, but 
they need to be minor enough to avoid strong morphologi- 
cal disruption. It is also possible that these galaxies are just 
approaching the cluster for the first time, and will eventually 
stop forming stars due to gas removal by the cluster envi- 
ronment. This would take them to the red sequence without 
severely disrupting their morphologies. 

• We observe no significant evolution of the intrin- 
sic colour scatter to z ~ 0.8 from the EDisCS clusters 
alone. This resu lt is consistent with previous studies (e.g. 
lEllis et al. Ill997l ). After expanding our sample with higher 
redshift clusters from the literature, we have still found no 
significant evolution in <7i nt up to z ~ 1.5. Moreover, in the 
wide range of cluster velocity dispersion (mass) of our sam- 
ple (100 < a v < 1300) the scatter does not seem to show 
any trend. Because our sample is strictly morphologically- 
selected, this implies that by the time cluster elliptical and 
SO galaxies achieve their morphology, the vast majority have 
already joined the red sequence. The only exception seems to 
be the very small fraction (< 7%) of faint blue early-ty pes. 

• Following the work of iBower. Lucev fe Ellis! J1992J '), we 
used the colour scatter to estimate the galaxies' formation 



time £f, defined as the time elapsed since the major episode 
of star formation. This allowed us to calculate the formation 
redshift zp for the early-type galaxy population in each clus- 
ter. Yet again, we measured no significant dependency of z-p 
on the cluster velocity dispersion. However, we found that 
zp increases weakly with cluster redshift within the EDisCS 
sample. This tren d becomes v e ry cle ar when the higher red- 
shift clusters from lMei et al. I (120091 ) are included. This im- 
plies that, at any given redshift, to have a population of 
fully-formed ellipticals and SOs they must have formed most 
of their stars ~ 2-4 Gyr prior to observation. That does not 
mean that all early-type galaxies in all clusters formed at 
these high redshifts. It means that the ones that we observe 
to already have early-type morphologies also have reason- 
ably old stellar populatio ns. This is partly a manifest ation 
of the "progenitor bias" (|van Dokkum fc FranxHl996T ). but 
also a consequence of the vast majority of the early-type 
galaxies in clusters (in particular the massive ones) being 
already red (i.e., already having old stellar populations) by 
the time the achieved their morphology. 

• Elliptical and SO galaxies show very similar colour 
scatter, implying that they have similar stellar population 
ages. If we assume that their observed properties are rep- 
resentative of the early-type cluster galaxy population at 
these redshifts, the scarcity of blue SOs indicates that, if 
they are t he de scendan ts of spirals whose star formation 
has c e ased (|Ara gon-Sala manca et al. 1120061 : iBedregal et al. I 
120061 ; iBarr et al. I [2007), the galaxies were already red 
when they became SOs, i.e. the parent spiral galaxies be- 
came red before loosing their spiral arms. T he red spi- 
rals f ound preferentially in dense environm ents (Wolf et al. I 
l2009l ; lBamford et al. Il2009l : iMasters et al. II2OO9T ) are the ob- 
vious candidate progenitors of these SOs. 

• Dividing the sample in two halfs by luminosity (or 
stellar mass), we find that the formation redshift zp (de- 
rived from the CMR scatter in each sample) is smaller for 
fainter galaxies than for brighter ones. This indicates that 
fainter early-type galaxies finished forming their stars later. 
Our results are also consistent with the observation that 
the cluster red sequence built over time with the bright- 
est galaxies reaching the sequence earlier than fain ter ones 
l|De Lucia et al. Il2004bl . I2007J : iRudnick et ai~ll2009h . 

• The CMR scatter analysis cannot constrain both the 
formation time tp and formation interval At simultaneously. 
However, its combination with the observed evolution of the 
CMR zero point, enabled us to conclude that the early-type 
cluster galaxy population must have had their star forma- 
tion truncated/stopped over an extended period At > 1 Gyr. 
Hence, the cessation of star formation was not synchronized 
for all the cluster early-type galaxies. 
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